Scattering of cosmic strings by black holes: loop formation 
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We study the deformation of a long cosmic string by a nearby rotating black hole. We examine 
- - - whether the deformation of a cosmic string, induced by the gravitational field of a Kerr black hole, 

may lead to the formation of a loop of cosmic string. The segment of the string which enters the 
, ergosphere of a rotating black hole gets deformed and, if it is sufficiently twisted, it can self-intersect 

^ ■ chopping off a loop of cosmic string. We find that the formation of a loop, via this mechanism, 

is a rare event. It will only arise in a small region of the collision phase space, which depends on 
^ |. the string velocity, the impact parameter and the black hole angular momentum. We conclude that 

D ■ generically, the cosmic string is simply scattered or captured by the rotating black hole. 
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PACS numbers: 04.20.-q, 04.70.Bw, 11.27.+d 



I. INTRODUCTION 
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■ Cosmic strings are one-dimensional topological defects [l|, y{ which could have been formed at the end of phase 
transitions associated with spontaneously broken symmetries in the early universe, via the Kibble mechanism |J|. 
Even though cosmic strings do not play a dominant role in structure formation, these objects are expected to be 
generically formed in the framework of supcrsymmetric grand unified theories [3- Cosmic microwave background 
measurements allow a small, but non negligible, contribution of strings 0, which can be realized in supcrsymmetric 
inflationary models 0. Recently, cosmic strings gained a lot of interest since it was pointed out that fundamental 
strings, Dl-branes and D(p — 2)-branes wrapping around {p — 3)-cycles in the compactified dimensions, could also 



I ' 

^ ' play a cosmic string-like role in cosmology 
rr| Cosmic string interactions can lead to a network of string loops. Intercommutations of two long strings in two 
^ I points, or self-intercommutations of a long string in one point, may chop off a loop of cosmic string. In what follows, 
' we explore another mechanism which may lead to the formation of cosmic string loops. More precisely, we examine 
, the possibility that the scattering of a long cosmic string by a rotating black hole can cause a deformation of the long 
d ' string, which may lead to the formation of a loop of cosmic string. As a segment of a long cosmic string enters the 
ergosphere of a rotating black hole, with rotation axis not parallel to the axis of the cosmic string, this segment may 
be deformed and, in some cases, it may get twisted. If the long string self-intersects, as a result of the deformation, 
a loop of cosmic string may be formed. 

The formation of loops of cosmic strings is important for the astrophysical implications of a string network. In 
particular, cosmic strings can be detected through the gravitational wave background produced by oscillating string 
loops H \ 

The scattering of a long cosmic string by a rotating black hole has been previously studied in the case of a black 
hole whose axis of rotation is parallel to the axis of the long string, concluding that the long string will be either 
scattered, or captured by the black hole. Which of the two regimes will take place depends on the impact parameter. 
The critical value of the impact parameter which separates the two regimes is g iven in Ref. [m. The effect of the 
scattering of a long cosmic string by a rotating black hole was studied in Ref. |l2l | , where it was shown that the string 
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is displaced in the direction perpendicular to its velocity by an amount which depends on the impact parameter |l2l| . 
In this article we generalize the study of Ref. |l2l | for different configurations of a long cosmic string approaching a 
rotating black hole. 

We conclude that the formation of a loop from the scattering of a long string by a rotating black hole is a rare 
event. It will only arise in a small region of the collision phase space, which depends on the string velocity, the impact 
parameter and the black hole angular momentum. More frequently, the cosmic string is either simply scattered or 
captured by the rotating black hole. 

In what follows we assume that the string thickness is negligible as compared to its length. Moreover, the string 
linear mass density is very small. Then, assuming that the string size is much larger than the Schwarzschild radius 
of the black hole, while the mass of the string is much smaller than that of the black hole, we can treat the cosmic 
string as a test object, and we can neglect its gravitational back reaction. 

The string is initially sufficiently far from the black hole, so it can be taken as a straight string. As it approaches 
the black hole, it will move in a Kerr geometry, however since the string is taken to be sufficiently long, the end 
regions of the infinitely long string will not be influenced by the presence of the gravitational field of the Kerr black 
hole. 

Let us consider the axis of rotation of the Kerr black hole to be along the direction A", in a system of coordinates 
(X,y, Z). A cosmic string being initially far away from the black hole is moving towards it. One can investigate the 
following setups: 

• The string is at first parallel to the axis of rotation of the Kerr black hole, and it moves towards the black hole 
with a velocity perpendicular to the axis of rotation. Thus, the cosmic string is originally along the A'-axis and 
its velocity is along the Z-axis. [We call this case (i).] 

• The string is at first perpendicular to the axis of rotation of the Kerr black hole and it moves towards the black 
hole with velocity parallel to axis of rotation. Thus, the cosmic string is originally along the Z-axis and its 
velocity is along the A'-axis. [We call this case (ii).] 

• The string is at first perpendicular to the axis of rotation of the Kerr black hole and it moves towards the black 
hole with velocity which is also perpendicular to the axis of rotation. Thus, the cosmic string is originally along 
the Z-axis and it moves towards the Kerr black hole with velocity along the J^-axis. [We call this case (iii).] 

Case (i) has been studied in Ref. what follows we will focus on cases (ii) and (iii). Case (iii) is the most 

promising configuration to lead to the formation of a cosmic string loop. 

We organize the rest of the paper as follows: In Section II we briefly discuss the string equations of motion in a 
curved background. We give the string equations of motion for a string in a general background geometry and then 
we derive them for a straight string far away from a black hole. We proceed with first the Newtonian and then the 
Lense-Thirring scattering of a straight string in the linearized approximation. In Section III we study the scattering 
of a long straight string by a Kerr black hole. Our numerical approach is presented in Section IV. We round up with 
our conclusions in Section V. 

We will use units where G = c = 1. 



The world history of a cosmic string can be expressed by a two-dimensional surface in the four-dimensional space- 
time, which is called the string worldsheet: 



where the worldsheet coordinates C^jC^ a-re arbitrary parameters chosen so that C*^ is time-like (= t) and is 
space-like (= cr). 

The string equations of motion, in the limit of a zero thickness string, are derived from the Nambu-Goto effective 
action which, up to an overall factor, corresponds to the surface area swept out by the string in space-time: 



II. STRING EQUATIONS OF MOTION IN A CURVED BACKGROUND 



A. String equations of motion 



(II.l) 




(11.2) 
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where 7 is the determinant of the two-dimensional worldsheet metric ^ab'- 
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7 = det(7ah) = ^e'"'e'"^ -fabled , lab = 9 iiuX^a^^^i, , 



(II.3) 



with g^^i, the four-dimensional metric. 

One can derive the same string equations of motion by using Polyakov's form of the action [l^ 



where hab is the internal metric with determinant h. 

Varying Eq. (|II.4p with respect to x^{(°-) we obtain the string equations of motion 



lab - ^habh"'^!, 



cd 



where 



(dynamical) 
(constraint) 



(11.4) 



(II.5) 



□ = 



-h 



(11.6) 



and r(^g. is the four-dimensional Christoffel symbol. 

We choose the gauge in which hab is conformal to the flat two dimensional metric r]ab ~ diag(— 1, 1). The dynamical 
and constraint equations, Eqs. (|II.5p . become 



701 = 5M''^!o^ri 
7oo + 111 = dt^i^x'^o^'^Q + x'^^x'^^ 






, 



(II.7) 



where □ = -9^ + (note that Oa = d/dA). 



B. Straight string moving in a weak gravitational field 

As long as a cosmic string is far away from a black hole, it just feels a weak gravitational field and thus the linearized 
analysis is sufficient. We present below the linear approach, which can be employed to furnish the initial conditions 
allowing us to perform a full numerical simulation for the case of a string being close to a rotating black hole, and 
thus limiting the size of the numerical simulation. This approach has been already used in Ref . j^^l ■ 

Assuming at first that there is no external gravitational field, the 4-dimensional space-time metric g^i/ is just the 
Minkowski metric r}^^. In coordinates {T^X^y,Z) and choosing signature (—,+,+,-!-) one can show that a solution 
of the equations of motion Ea. ljII.5|) read 

x^ = Xq (t, cr) = (r cosh(/3), r sinh(/3) + Ao, cr) 
or x^ = Xq (t, cr) = (rcosh(/3), 5, Tsinh(/3) +3^0,0") , 

hab = ??afc = diag(-l, I) . (II. 8) 

Equation Ijll.Sb ,') refers to a long straight string along the Z-axis, being initially at a::^(0, a) = (0, Aq, 6, cr), and moving 
with velocity v = tanh(/3) along the direction X. We consider 6 > 0, Aq < 0. This configuration corresponds to case 
(ii). Similarly, Eq. IjlLSb ') refers to a long straight string along the Z-axis, being initially at 2:^(0, cr) = (0, 6,3^07 f), 
and moving with velocity v = tanh(/3) along the direction y. We again consider 6 > 0,3^o < 0. This configuration 
corresponds to case (iii). 

As a second step we consider the case of a weak gravitational field 

in which moves a straight long string 



(II.IO) 
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Sfiiy,x^{^) denote a small metric perturbation, and a small perturbation along a straight cosmic string, respectively. 
In first order in 0^1, and x'^ , the equations of motion, Ea. (|11.5|l . read 



ax>' + T'^p{xo)xiy^,7j'''' - (dynamical) 

dxn dx" dxf^ dxK ^ dxii dxn ^ . . , 

Vt^y-^-^ 'riif.y-^ ^ VOt,„-^^T— = (constramt) 

OT aa OT oa or oa 



/dx^dx" dx'idx''\ ^ /dx'r!.dx^ dx[; dx^ ^ , , , , 

{^^ + + -7^^ + =0 constraint . 11.11 

\ OT OT Oa OCT J \ OT OT Oa oa ' 



The above equations describe the motion of a long straight cosmic string, being located far from a Kerr black hole. 

In the weak field approximation, the gravitational field produced by a black hole of mass M and angular momentum 
J, rotating around the A'-axis, is 

ds^ = - (^1 - dT^ +(^1 + ^^ {dx^ + dy^ + dz^) - ^{ydz zdy)dT, (11.12) 

where TZ"^ = X'^ + y^ + Z^ . 

One can therefore write the small perturbation on the gravitational field, in terms of the Newtonian and the 
Lense-Thirring 0,0] parts, as 

4 / 

= 2^5^, + —5l^e,)^^,X'- , (11.13) 

where (p is defined as (p ~ M/TZ and eap-ys stands for the antisymmetric tensor. The second term of the r.h.s of 
Eq. pi.l3b l is due to the rotating black hole around the A'-axis. 
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C. Newtonian scattering of a straight string in the Hnearized approximation 

Let us first consider the Newtonian scattering of a straight cosmic string. We first examine case (ii). The dynamical 
equations read 

x° = -2sinh(/3)cosh(/3)(^,i 
x^ = Q 

x? = -2sinh^(/?)v3,2 

X? = -2cosh2(/3)(^,3 , (11.14) 
and the constraint equations are 

irjV - cosh(/3)x°^ + sinh(/3)x^^ = 

ij^^ - cosh(/3)x", + sinh(/3).i'^^ = -2v9Cosh2(/3) ; (11.15) 

we use the notation >V,a = dW jdA. 

The dynamical equations imply that for a given time the cosmic string lies in the (3^, Z)-plane. The string deflection 
T) along the plane orthogonal to the velocity of the string, given in Ref. i is 

V = b- 27rA/sinh(/3) . (11.16) 

It implies that while at later times segments of the cosmic string which are far away from the rotating black hole keep 
moving in the (y = b, Z)-plane, since they are unaffected from the gravitational field of the black hole, the central 
part of the string being by definition closer to the location of the black hole will be moving in the (3^ = V, Z)-plane 
instead. 
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We then study case (iii). The dynamical equations read 



and the constraint equations are 
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9x2 ~ 


^ 9^2 
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9t2 ' 


^ 9cr2 



x° = -2sinh(/3)cosh(/3)(/3,i 



= -2sinh^(/3)v3,i 


= -2cosh2(/3)(^,3 , (11.17) 



xjV - cosh(/?)x°^ + sinh(/3);x^^ = 

- cosh(/3)x°, + sinh(/3).i:;V = -2ipcosh^{p) . (11.18) 



D. Lense-Thirring scattering of a straight string in the linearized approximation 



Let us now proceed with the Lense-Thirring scattering of a long straight string, moving in the gravitational field of 
a black hole of mass M and angular momentum J, in the linearized approximation. 
Starting with case (ii), the dynamical equations read 



9^ 
"9t2 

9^ 
'9t2 

9^ 
'9t2 

9^ 

9t2 



9^2 
&^ 

92 \ 
9a2 ) 

&^ 
9ct2 



QJba 



7^5 



= 



6J 

— ^(Tsinh(/3) + a;o)crcosli(/3) sinh(/3) 



6J 



(Tsinh(/3) + a;o)&cosh(/3) sinh(/3) 



and the constraint equations are 



x\ - cosh(/3)x"^ + s\nh{(3)x\ = 



2J, 



where 



Choosing initial conditions 



- cosh(/3)x°^ + sinh(/3)x^^ + irj'^ = ^focosh(/3) 
n = v/(t sinh(/3) + XoY + 62 + ^2 . 



9a;2 



a; r=o 



|t=0 







9t 
9r 



(n.i9) 



(n.20) 



(11.21) 



(n.22) 
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the dynamical equations can be solved analytically. The solution reads 
2Jb 



sinh^(/3)T + Xo sinh(/3) - a sinh^ {(3)t + Xq sinh(/3) + a 



Air + a)n Air - a)n 

T + (J — Xq sinh(/3) t — a — Xq sinh(/3) " 



x-" = 2J6cosli(/3) 



A{t + (7)S{t + a) A{T-a)S{T-a) 

sinh(/3)[sinh(/3){62 + + err} + aXo] sinh(/3)[sinh(/3){62 + - ar] - aX^] 



A{t + (T)7^ A{t - a)Tl 

sinh(/3)[sinh(/3)62 + sinh(/3)(r + af + A'o(t + cr)] sinh(/3)[sinh(/3)&2 + sinh(/3)(r - af + Xq{t - a)\ 



X? = 2J6cosh(/3) 



A{T + a)S(T + a) A(T-G)S(T~a) 
sinh^ (/3)t + Xq sinh(/3) + a sinh^ (/3)t + Xq sinh(/3) - cr t - cr - A'o sinh(/3) 
A(t - a)n ^ ;4(7T^0^^ ^ ^(r - (7)5(t - a) 

T + a ^ Xq sinh(/3) r + cr t — cr 



A{T + a)S{T + a) iXg+y^)S{T + a) [X^ + y^)S{T - a) 

where 



7e = Jsmh^{p)T^+2smh{l3)TXQ+Xg+y^ + cT^ 



= 



and the constraint equations arc 



where 



(11.23) 



siu) = ^xi + y^ 

A{u) = yl[l + sinh(/?)2] + [Xq + sinh(/3)w]2 . (11.24) 
Let us proceed with case (iii). The dynamical equations read 

= -^(rsinh(/3)+2/o)cosh2(/3) 

+ = ^ cosh(/3) smh(/3) 



2 J 

x^^ - cosh(/?)S°^ + sinh(/?)s^^ = -—cr cosh(/3) sinh(/?) 

2 7 

cosh(/3)x°^ + sinh(/3)xf^ + xjV = ^6cosh(/3) , (11.26) 



7^ = V&2 + (^sinh(/3)+:^^o)2+a2 . (11.27) 
The solution is far more complicated than the one for case (ii)^. 

III. SCATTERING OF A STRAIGHT STRING BY A KERR BLACK HOLE 

We now turn to the full non-linear analysis. The Kerr metric, can be presented in the quasi-Cartesian coordinate 
system (T, X, y, Z) by using the following shift from the Boycr-Lindquist system (t, r, 9, (p): 

r = n + M, (III.28) 



^ Instead of writing down explicitly the quite tedious solution, wc give in an appendix the method we follow to solve the system of 
equations. 
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and 

r = t 

y = 7?. sin cos (/) 

Z ^ 7^sin6lsin<^!) , (III.29) 
In these coordinates, the rotation axis of the black hole is the X axis and the metric takes the form 



(3^2 + 22)7^2A2 ' 7^2 7e2(3;2 + 2:2) 

A2 /A'2 J^2^22^ ,^^2 ,2^2/1 1 



7^2 ' 7^2 + + (III.30) 

where M, J denote the mass and angular momentum, respectively, of the black hole, withe J = aM for |a| <M; A 
and p are defined as 

^2^2 

= (7^ + A/)^ + —5- and = (7^ + M)^ + - 2M (Tl + Af) . (III.31) 
Far from the origin of the coordinate system, the metric takes the form given in Eq. (|II.12|) . 



IV. NUMERICAL APPROACH 



We use the quasi- Cartesian coordinate system given previously. To study the case (ii), or (iii), we let evolve a 
string initially parallel to Z axis with initial velocity along the X axis (case (ii)), or the y axis (case (iii)). We focus 
on the more promising case of a maximally rotating black hole (a = M) and we set the mass equal to unity; we can 
recover other values by rescaling the coordinates. 

The Kerr metric is given analytically in Eq. (|III.30|I and the Christoffel symbols are obtained by taking the deriva- 
tives of the metric. In our simulations we compute the derivatives as a finite difference for a small variation of the 
coordinates. We discretise a piece of world-sheet with steps Ac and At. The string position is obtained as the value 
of the fields X^ia^ r) on the world-sheet. They are computed using Eq. Ijll.7|) and performing a numerical integration 
along the r direction. We need initial and boundary conditions: the initial conditions arc the value and t derivatives 
of X^[a^ 0), V^ip-: 0) = drX^{a, 0) on the slice r = 0. We compute the boundary conditions as the free evolution of 
the end of the world-sheet'^. 

In our simulations the initial conditions are the values of the X^((t, 0), ^^(cr, 0) at the point ai, i G {0; ...;4000} 
given by 



for the case (ii), and 



X^(cr,,0) = {Q,-Xo,b,{i- 2000)6Z) , 
V^{a,,0) = (cosh(/?),sinh(/3),0,0) , 



X'^ia^O) = (0,6,-3^o,(i- 2000)(5Z) , 
F^(a,,0) = (cosh(/3),0,sinh(/3),0) , 



(IV.32) 



(IV.33) 



for the case (iii). In the above, b is the impact parameter, 6Z = 0.1/A/, f3 denotes the velocity parameter and Xq (J^o) 
is the distance between the cosmic string and the black hole at the beginning of the simulation. In order to avoid 
numerical instability we set the proper time step At = 10^'^. 

For a cosmic string which is deformed by a Kerr black hole but with insufficient twist, so that a loop is not formed 
when the string comes close to the black hole, we stop the simulation when the string starts drifting away from the 



In a more accurate simulation both initial and boundary conditions should be computed using the weak field evolution. 
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black hole. The reason for this is that the deformation of the central segments of the cosmic string will be propagated 
towards the end parts of the string, leaving behind a wiggly, but on the large-scale straight, cosmic string. We also 
stop the simulation if we encounter an infinity, denoting that a point along the string crosses the horizon of the black 
hole. 

To display the string position we produce equal-time slices. Such slices are obtained by searching all the world-sheet 
points (cr, t) such that X°((T, t) = T for fixed T and plotting the {X, Y, Z). 



200- 
100. 




40 X 



Y 

FIG. 1: Timeslices for the case (ii) with /3 = 0.5 and 6 = 5. We show the slices for time T = 25, 37.5, . . . , 100. 




FIG. 2: Same as Fig. but with the X motion subtracted in order to show the string deformation. 



V. CONCLUSIONS 

We have studied the deformation of a long cosmic string entering the ergosphere of a Kerr black hole. The aim of 
this study is to investigate how often such a deformation can be accompanied by a sufficient twist, such that there is 
a string self-intersection leading to a loop formation. 

Among the configurations we have studied, the case of a long cosmic string being at first perpendicular to the axis 
of rotation of the Kerr black hole and moving towards it with velocity perpendicular to the rotation axis, is the most 
prominent to lead to loop formation. [It is our case (iii).] 

A numerical investigation has shown us that loop formation is a rather rare event, while in general the cosmic string 
will be just scattered from the black hole or captured by it. 
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100. 
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FIG. 3: Timeslices for the case (iii) with (3 = 0.5 and b — 3.3. We show the shoes for time T = 44, 55, . . . , 110. 



200-1 




X Y 



FIG. 4: Same as Fig. |3]but with the y motion subtracted in order to show the string deformation. Note that the center of the 
string has been twisted of about 90 degrees. 
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Appendix 



Rather than giving explicitely the lengthy formulae for the solutions of the dynamical equations, Eq. Ijn.25|l . of case 



(iii) with constraint equations Eq. (|II.26p . we describe here the resolution method of such equations. 
The form of the equations is 

9^<^(a,r)-a,V(^,r)=/(a,r) . 
The corresponding homogeneous equation is 



(V.34) 
(V.35) 



The general solution of Eq. (|V.34|I is the sum of the general solution (ph of the homogeneous equation Eq. ljV.35|l 
and a particular solution tppart of Eq. (|V.34p . 
The general solution tph is just 



tphlo", t) = a{a + t) + b{(T — t) , 



wih arbitrary functions of one variable a and b. 



It is possible to write a particular solution (/'part a-s a double integral of the right hand side of ljV.34|) 



<Ppart = / .f{x,y) dx dy , 
I A 



(V.36) 



(V.37) 



where the domain of integration A is a triangle shown in the figure below. 





y 


(J + r 












a — r 





X 



The particular solution (ppart may be written in more than one way. Two are useful 



rr-cr+y ^ 

dy / dxf{x,y) + - 

Jo ^ 



(T+r r-T+a-y 

dy / f{x,y) 
Jo 



(V.38) 



and 
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2 /•T /■a+T—x 

<<5part(o-,r) = - / dx I dyf{x,y) (V.39) 

Expression (|V.38|I [rcsp. (|V.39|l ] is useful when f{x,y) may be written simply as a derivative dxg{x,y), [rcsp. 
dy g{x,y)]. Only one integration is then to be performed. 

Of course equation (|V.34|) is complemented with boundary conditions which will determine the unknown functions 
a and b. 



